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ABSTRACT 

Although it is well recognized that gamma-ray burst (GRB) afterglows are obscured and reddened by dust in 
their host galaxies, the wavelength-dependence and quantity of dust extinction are still poorly known. Current 
studies on this mostly rely on fitting the afterglow spectral energy distributions (SEDs) with template extinction 
models. The inferred extinction (both quantity and wavelength-dependence) and dust-to-gas ratios are often in 
disagreement with that obtained from dust depletion and X-ray spectroscopy studies. We argue that this discrep- 
ancy could result from the prior assumption of a template extinction law. We propose an analytical formula to 
approximate the GRB host extinction law. With the template extinction laws self-contained, and the capability 
of revealing extinction laws differing from the conventional ones, it is shown that this is a powerful approach in 
modeling the afterglow SEDs to derive GRB host extinction. 

Subject headings: dust, extinction — gamma rays: bursts 



1. INTRODUCTION 

In addition to the Galactic foreground extinction, GRBs and 
their afterglows are subject to extinction caused by the dust 
within their host galaxies. Evidence for this includes — 

• "Dark bursts" - an appreciable fraction of GRBs with 
X-ray and/or radio afterglows lack an optical afterglow 
(Jakobsson et al. 2004)0 A natural explanation for 
dark bursts is that they lie behind significant obscur- 
ing dust columns in their host galaxies which effectively 
suppresses the optical light [although some dark bursts 
may be intrinsically faint or occur at high redshifts (say, 
z > 5) where the Lya break has moved through the op- 
tical bands, leading to absorption of the optical light by 
the Lya forest]. Indeed, Schady et al. (2007) found that 
the X-ray afterglows of GRBs not detected by UVOT 
were more affected by extinction than those of GRBs 
with detected UVOT counterparts. The recent detec- 
tion of the near infrared (IR) afterglows of some GRBs 
(which would have been considered as "dark bursts" 
since their afterglows were not detected in any bluer 
bands) provides another piece of evidence for dust ob- 
scuration (e.g. see Jaunsen et al. 2008, Tanvir et al. 
2008). 

• Reddening - some GRB afterglows with low redshifts 
appear very red, due to effects of extinction - ultravio- 
let (UV)/visible light is extinguished more by dust than 
red light (e.g. see Klose et al. 2000, Levan et al. 2006). 
Dust reddening is also indicated by the significant de- 
viation of the optical/near-IR spectral energy distribu- 
tions (SEDs) of many afterglows from that expected 
from standard models. Also because of dust reddening, 
the Balmer line ratios in the spectra of some GRB host 
galaxies (e.g. see Djorgovski et al. 1998), known as the 



the standard Case B recombination, which are fairly in- 
dependent of physical conditions (Osterbrock & Ferland 
2006). 

• Depletion - dust-forming heavy elements such as Si and 
Fe were found to be substantially depleted from the gas 
phase in some host galaxies (e.g. see Savaglio et al. 
2003). This indirectly shows the presence of dust in 
GRB host galaxies since the missing heavy elements 
must have been locked up in dust grains. 

• Connection between long GRBs and massive stars - 
there are multiple strong lines of evidence that long- 
duration (> 2 s) GRBs are associated with the death of 
massive stars, occurring in regions of active star for- 
mation embedded in dense clouds of dust and gas (see 
Woosley & Bloom 2006). 

A precise knowledge of the extinction (quantity, wavelength- 
dependence) and the nature (size, composition, and quantity) of 
the dust in GRB host galaxies is crucial for 

• Correcting for the extinction of afterglows from X-ray 
to near-IR wavelengths to derive their intrinsic lumi- 
nosities - this is particularly important for studying the 
luminosity distribution of GRB afterglows and their in- 
trinsic SEDs (e.g. see Kann et al. 2008); 

• Constraining the nature of the GRB progenitors (i.e. 
collapsing massive stars or merging neutron stars) - if 
long-duration GRBs are indeed linked to the collapse of 
massive stars, it is most likely that their optical and near- 
IR afterglows will suffer from significant attenuation in 
the star-forming molecular clouds heavily enshrouded 
by dust - the birth place of these short-lived (~ 10 6 yrs) 
massive stars; 

• Tracing the physical conditions of (and processes oc- 
curring in) the environments where GRBs occur which 



Balmer decrement, deviate from the expected ratios for 
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hold clues for understanding the mechanism for making 
a burst, e.g., a fiat or gray extinction law for GRB host 
galaxies would imply a dense circumburst environment 
where dust undergoes coagulational-growth or a prefer- 
ential destruction of small grains; and 

• Probing the interstellar medium (ISM) of high-redshift 
galaxies and the cosmic star formation history - because 
of their intense luminosity which allows their detection 
at cosmological distances, GRBs are a powerful tool to 
study the star formation history up to very high red- 
shifts; e.g., the dust and extinction properties of GRB 
hosts would help understand the nature of dark bursts 
and the dark burst fraction which would place important 
constraints on the fraction of obscured star formation in 
the universe (e.g. see Djorgovski et al. 2001, Ramirez- 
Ruiz et al. 2002). 

However, our current understanding of the dust extinction in 
GRB host galaxies is still very poor. Existing studies on this 
often draw conclusions in conflict with each other (see $2] for 
details). We argue that this could be caused by the prior adop- 
tion of a template extinction law in fitting the observed GRB 
afterglow spectra to derive dust extinction We propose 

in this work an alternative, robust method based on an analyti- 
cal formula which can restore the widely adopted template ex- 
tinction laws ($3). For illustration, we apply this approach to 
GRB000301C and GRB 021004 (©. We demonstrate in g5] 
the uniqueness of the derived extinction laws. The robustness 
of this approach will be discussed in a separate paper (Liang & 
Li 2008a) in which the afterglow SEDs of > 50 GRBs of a wide 
range of properties are successfully modeled and for which the 
inferred extinction curves are diverse, with some differing sub- 
stantially from any of the template extinction curves. 

2. CURRENT STATUS 

At present, the amount of extinction (usually the rest-frame 
visual extinction Ay r ) and the wavelength-dependence of the 
extinction ("extinction curve" or "extinction law"; A\/Ay or 
A„/Ay if expressed in frequency) are commonly derived by 
fitting the UV, optical, and near-IR afterglow photometry (F„; 
with the Galactic extinction corrected) with a power-law model 
(oc v~P\ approximating their intrinsic spectra) reddened by an 
assumed, template extinction law A v jAy 



F„ = F (v/Hz) ^exp 



1.086 A 



v, 



(1) 



where (3 is the intrinsic power-law slope of the afterglow, F a is 
a normalization constant (normalized to the overall afterglow 
flux level), A(i +Z )i, is the rest-frame extinction, and z is the GRB 
redshift. The factor of "1.086" in eq.(fl]i arises from the conver- 
sion of extinction (in magnitude) to optical depth. As a priori, 
six template extinction laws have been widely adopted in the 
literature to derive the dust extinction of GRB hosts: (1) a sim- 
ple power-law Ax /Ay ~ A~ 7 or even just a linear function of 
inverse wavelength Ax/Ay ~ A" 1 ("Linear" thereafter); (2) the 
Milky Way (MW) extinction curve (with a prominent bump at 



2175 ) characterized by Ry, the total-to-selective extinction ra- 
tio (the Galactic average value is Ry ~ 3.1); (3) the feature- 
less Small Magellanic Cloud (SMC) extinction curve which 
steeply rises with inverse wavelength from near-IR to far-UV 
(Ax/ 'A v ~ A" 1 2 ); (4) the Large Magellanic Cloud (LMC) curve 
being intermediate between that of the MW and the SMC; (5) 
the featureless "Calzetti" attenuation law for the dust in local 
starburst galaxies (Calzetti et al. 1994)0 and (6) the relatively 
flat "Maiolino" extinction law for the dust in the dense circum- 
nuclear region of AGNs (Maiolino et al. 2001) where the dust 
size distribution is skewed toward large grains (see Fig. 1). 

To our knowledge, exceptions to the "template" extinction 
approach described here are that of Chen et al. (2006) and Li 
et al. (2008), both of which were based on the fireball model. 
The latter approach is limited to bursts of which the X-ray and 
optical decay indices are the same. In most studies (which as- 
sume an extinction template) a SMC-type extinction curve is 
preferred. This is probably because the 2175 extinction fea- 
ture (which is prominent in the MW and LMC curves) is rarely 
seen in the afterglow spectra of GRBs. So far, its possible de- 
tection is only reported in four bursts: GRB 970508 (Stratta 
et al. 2004), GRB 991216 (Kann et al. 2006, Vreeswijk et al. 
2006), GRB 050802 (Schady et al. 2007), and more definitely 
GRB 070802 (Kruhler et al. 2008; A. Etfasdottir et al. 2008, in 
preparation). 

However, some studies favour a much flatter or even gray ex- 
tinction curve (e.g. see Savaglio et al. 2003, Savaglio & Fall 
2004, Stratta et al. 2005, Chen et al. 2006, Li et al. 2008, Per- 
ley et al. 2008). With a SMC-type curve, the amount of vi- 
sual extinction Ay or reddening] derived by fitting the after- 
glow photometry tends to be small since the SMC curve rises 
so rapidly with A -1 that a small Ay would imply a large UV 
extinction. This may explain the finding of "a strong clustering 
toward low extinction (Ay < 0.2 mag)" in a detailed study of 19 
GRBs by Kann et al. (2006), and later by Kann et al. (2008) 
for 15 GRBs. In contrast, for a flatter extinction law like that of 
Calzetti, Maiolino, MW with R v > 4, or that derived by Chen et 
al. (2006), Li et al. (2008) and Perley et al. (2008), a relatively 
large Ay is often obtained. 

The visual extinction Ay can also be inferred from the dust 
depletion method based on the gas-phase heavy-element abun- 
dances estimated from the afterglow optical absorption spec- 
troscopy (Savaglio et al. 2003; Savaglio & Fall 2004). This 
analysis assumes both the dust depletion pattern and the visual 
extinction per unit dust column Ay /A/dust of GRB hosts to be the 
same as that of the MW. It is quite possible that GRB hosts may 
have a different depletion pattern and/or a different Ay /AWst 
conversion factor. The latter could result from a dust composi- 
tion or size distribution differing from that of the MW. 

One can also derive Ay from the neutral hydrogen col- 
umn density A^h derived from Lya absorption (Hjorth et al. 
2003) or the equivalent A^h obtained from soft X-ray absorp- 
tion (mostly from oxygen K-shell absorption; Galama & Wi- 
jers 2001, Stratta et al. 2004, Watson et al. 2006). There is a 
puzzling discrepancy between the optical reddening E(B-V) 
derived from the afterglow SED fitting and the visual extinc- 
tion Ay inferred from the dust depletion analysis or from A^h 



6 We should note that recent Spitzer observations in the near- and mid-IR argue against GRB hosts being strongly starbursting galaxies (Le Floc'h et al. 2006), 
although their morphological and average radio/submillimeter properties suggest that they are likely massive and actively star-forming galaxies (Berger et al. 2003; 
Conselice et al. 2005). 

7 Reddening is usually expressed as E(B-V) = Ag-Ay = Ay/Rv, where Ag is the extinction at the B band (\g 4400). By definition, gray dust (for which the 
extinction is just weakly dependent on A) is characterized by small reddening E(B— V) and large Ry. Apparently, for gray dust, a small reddening does not necessarily 
imply a small extinction since Rv can be large. 
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IR/visible extinction and the 2175 extinction bump. We call 
this the "Drude" approach since Eq.(f2| looks like a sum of 
Drude functions. As shown in Figure Q] this formula, with the 
free parameters cj (j = 1,...,4) adjusted using the Levenberg- 
Marquardt minimization algorithm (Press et al. 1992; see Table 
Q]), can reproduce the extinction curves widely adopted as tem- 
plate extinction laws in GRB afterglow SED modeling, clearly 
demonstrating the advantages of the proposed formula over any 
template extinction laws with a fixed wavelength-dependence 
shape: with the widely-adopted conventional extinction laws 
self-contained in Eq.© and the capability of revealing extinc- 
tion laws differing from the conventional ones, the proposed 
formula is more flexible and more powerful in modeling the 
afterglow SEDs. Indeed, as shown in Liang & Li (2008a,b), 
dust reddening models based on this formula nicely reproduce 
the observed afterglow SEDs of distant GRBs at z > 4 (in- 
cluding GRB 050904 at z ~ 6.3) and that of the "troublesome" 
GRB 061 126 (Perley et al. 2008) without resorting to an exotic 
extinction law. 

4. GRB 000301C AND GRB 021004: TEST CASES 

We apply the above-described technique (|3]) to the optical 
afterglows of GRB000301C at z » 2.04 (Jensen et al. 2001) 
and GRB 021004 at z « 2.33 (Fynbo et al. 2005). They are 
selected mainly because they are among the best-observed in 
terms of sampling in the time domain, and multiwavelength 
coverage. We fit their broadband SEDs using Eqs. fl 1I21 > with 
(3, Ay, c\, C2, and c 3 allowed to vary as free parameters [F is 
not really a free parameter; for a given set of (/3, Ay, c\, C2, 
C3), F a is uniquely determined by the overall flux level. There- 
fore, in the SED modeling we fit five free parameters to the six 
(seven) data points of GRB 00030 1C (GRB 021 004)] We de- 
rive the best-fit parameters based on the Levenberg-Marquardt 
minimization algorithm (see Tabled. As shown in Figure|2]for 
GRB 00030 1C and in Figure|3]for GRB 021004, almost perfect 
fits to the observed SEDs are achieved through this approach. 
The inferred extinction curves differ substantially from any of 
the template extinction laws. 

5. DISCUSSION 

We have also fitted the afterglow SEDs of GRB 00030 1C and 
GRB 021004 in terms of the MW, SMC, LMC, Calzetti, and 
"linear" template extinction curves (see Tableland Figs. l2l31 >. 
Since for a given template extinction law the wavelength- 
dependence of the extinction A\/Ay is fixed, we are now left 
with only three parameters: F a , (3, and Ay. The models based 
on the MW and LMC extinction laws could not fit the observed 
SEDs at all. This is because the 2175 extinction feature which 
is prominent in the MW and LMC curves is absent in the SEDs 
of GRB 00030 1C and GRB 021004. In contrast, the SMC and 
"linear" models closely fit the afterglow SEDs of these two 
bursts, better than the Drude model proposed here as measured 
by X /Nd.o.f. (see Table©. 

While the Drude model has three more parameters than the 
SMC and "linear" models, the quality of the fitting of the Drude 

8 Reichart (2001) proposed a seven-parameter formula for the dust extinction curve Ax/Ay of GRB hosts based on the expressions of Cardelli et al. (1989; "CCM"; 
for A > 0.3 frni) and of Fitzpatrick & Massa (1990; "FM"; for 0.1 fim < A < 0.3 jum), The problem with the Reichart (2001) formula (see his eqs. 61,66) is that the 
CCM expression is only valid for the Galactic extinction curves, it is not suitable for the SMC or LMC extinction (Gordon et al. 2003). Therefore, if a GRB host 
happens to have a SMC- or LMC-type extinction law, models based on the Reichart (2001) formula will not be able to restore the true extinction. 

9 We set c 4 = based on a visual inspection of the observed SEDs which clearly suggest the absence of a 2175 feature (see Figs. |2|3> . With c 4 treated as a free, 
positive parameter, even the best fits (given by c 4 « 0.0034,0.0018 for GRB 000301C and GRB 021004, respectively; for comparison, c 4 « 0.051,0.039 for MW and 
LMC, respectively) are not as good as that provided by models with c 4 = 0. We place an upper limt of c 4 w 0.015 (0.0073) for GRB 000301C (GRB 021004). 



measured from the Lya or X-ray absorption spectra, with the 
former considerably smaller than the latter. 

This discrepancy problem could be alleviated if one invokes 
a flat or gray extinction law. This is because (1) gray dust 
(> 1 /im) characterized with an extinction curve weakly depen- 
dent on A in the optical/UV could produce a high Ay but little 
reddening (see Footnote-|7]l, and (2) per unit mass gray dust is 
not as effective as submicron-sized dust in absorbing and scat- 
tering optical light so that the Ay /Nfast conversion factor for 
gray dust is smaller than that of the MW dust (with a typical 
size of ~ 0.1 fim; see Li 2008). The latter would imply that 
the methods based on dust depletion (see above) and Lya/X- 
ray absorption may overestimate Ay if the dust size distribu- 
tion of GRB host galaxies is indeed biased toward large grains, 
as a result of dust coagulational growth in the dense circum- 
burst environments or preferential destruction of small dust by 
GRB emission (e.g. see Waxman & Draine 2000, Fruchter et 
al. 2001,Perna etal. 2003). 

However, we should stress that the gray extinction hypothe- 
sis should not be considered as the only solution to the discrep- 
ancy problem (after all, SMC-type or even steeper extinction 
laws were derived for the hosts of some GRBs; see Liang et al. 
2008a). Indeed, the mismatch between the X-ray-derived Ay 
and that derived from the optical SED modeling could be at- 
tributed to physically separate X-ray and optical emission re- 
gions (e.g. see Prochaska et al. 2006, Watson et al. 2007). 
Prochaska et al. (2006) argued that for GRB 051 111, the X- 
ray opacity comes from dust-free gas that is very local to the 
GRB (^lpc), while they placed a lower limit of >50pc on 
the host galaxy absorption systems from the GRB. It has also 
been argued that Ay is probably probing the dust outside of the 
dense molecular cloud around the GRB, since all dust within 
the cloud is likely to have been obliterated by the burst (e.g. see 
Perna & Lazzati 2002, Prochaska et al. 2007). Moreover, if the 
dust depletion pattern of GRB hosts is different from that of the 
Milky Way, the discrepancy between the depletion-derived Ay 
and that from the optical SED modeling could be alleviated. 

3. OUR APPROACH 

In view of the shortcomings of the prior assumption of a tem- 
plate extinction law (see $2| and guided by Pei (1992), we pro- 
pose a simple formula containing four dimensionless parame- 
ters (ci, C2, C3, and c 4 ) for the wavelength-dependence of the 
extinction for the dust in GRB host galaxies, instead of adopt- 
ing any known extinction laws (see © as a template, 

c\ 



Ax/ Ay = 



(A/0.08) C2 + (0.08/A) C2 +c 3 
233 [l-c 1 /(6.88 C2 +0.145 C2 + c 3 )-c 4 /4.60] 
(A/0.046) 2 + (0.046/A) 2 + 90 

C4 

1.95 



(2) 



(A/0.2175) +(0.2175/A) 

where A is in /xm0 While the first term in the right-hand 
side of Eq.© represents the far-UV extinction rise, the sec- 
ond term and the third term respectively account for the near- 
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model is even not as good as that of the SMC or "linear" model. 
Then, why do not we simply adopt the SMC or "linear" model? 

First of all, we should note that there are no physical rea- 
sons for a prior assumption of a known extinction law, either 
that of the SMC, LMC, "linear" or MW: the composition and 
size distribution (and therefore the extinction law) of the dust in 
the dense circumburst clouds of GRB hosts with a wide range 
of metallicities and evolutionary stages are not expected to re- 
semble that of the MW, LMC, or SMC (e.g. see Dwek 2005). 
In literature, a SMC-type extinction is often assumed for low- 
metallicity environments. However, there is no physical basis 
for this (except the lack of grain growth in these regions because 
of the lack of raw dust materials - the SMC dust, on average, is 
substantially smaller than that of the Milky Way [see Weingart- 
ner & Draine 2001]). Moreover, it is known that the GRB hosts 
have a wide range of metallicities. Indeed, the reasons why the 
MW, LMC and SMC laws are often used for GRB afterglow 
SED modeling are mainly (1) little is known about the extinc- 
tion laws of other galaxies, and (2) the Pei (1992) formula for 
the MW, LMC and SMC extinction laws is numerically conve- 
nient for computer implementation. 

Second, although the SMC-type extinction is preferred in 
most of the present afterglow SED modeling studies, only the 
Drude approach is capable of reproducing the SEDs of those 
reddened by gray extinction or by non-conventional extinction. 
Indeed, it was shown that the afterglow SED of GRB 050904 
at a redshift of z ~ 6.3 cannot be explained by dust reddening 
with any of the conventional (MW, SMC, Calzetti) extinction 
curves; instead, it can be well reproduced by invoking the ex- 
tinction curve inferred for a distant quasar at z = 6.2 (Maiolino 
et al. 2004), suggesting that the properties of dust may evolve 



beyond z = 6 (Stratta et al. 2007). 

Third, the Drude model would at least complement the mod- 
els using template extinction curves, particularly for those 
bursts for which the Drude model gives a larger x 2 /Nd.o.{. (but 
still fits the observed SEDs well). Given that the derived ex- 
tinction Ay and the intrinsic spectral slope (3 differ appreciably 
among different approaches (see Tabled, the SMC model (and 
other models) should be used along side with the Drude model 
to gain insight into the "true" extinction and the "true" spectral 
slope. 

We finally demonstrate the uniqueness of the extinction curve 
inferred from the Drude approach. To this end, we generate 
three sets of afterglow "photometry data" by reddening the in- 
trinsic afterglow spectrum F v (piy) = 5.2 x 10 8 (V/Hz) ' of a 
burst at z ~ 2 respectively with three template extinction laws: 
MW, SMC, and Calzetti, each with A v = 0.5 mag. We then ap- 
ply the Drude approach to these three sets of artificially-created 
GRB afterglow data. As shown in Figure[4] we uniquely re- 
store the MW, SMC, and Calzetti extinction laws: the inferred 
extinction curves are almost identical to that used to redden the 
intrinsic spectrum (the derived parameters [see Table[3| are es- 
sentially the same as those tabulated in Table[TJ. 
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FIG. 1 . — Extinction laws widely adopted as "templates" in GRB host extinction studies: the SMC law (upper panel: dashed black line; Pei 1992), the LMC law 
(upper panel: dashed blue line; Pei 1992), the linear A\ oc A - ' law (upper panel: dashed cyan line), the MW Galactic average extinction law (Ry = 3.1; lower panel: 
dashed black line; Pei 1992), the MW extinction law with Ry = 4.5 for dense clouds (lower panel: solid red line), the Calzetti starburst attenuation law (lower panel: 
dashed blue line), and the Maiolino law for AGN dust tori (lower panel: dashed cyan line; just like that of the MW with Ry = 4.5). Also shown are the "Drude" fits 
to these "template" extinction laws: SMC (upper panel: solid green line), LMC (upper panel: solid red line), Linear (upper panel: solid magenta line), MW with 
Ry = 3.1 (lower panel: solid green line), and Calzetti (lower panel: solid magenta line). 



Table 1 

"Drude" fits to known extinction curves for A = 0.1-1 widely adopted as "templates" in modeling GRB afterglow SEDs to derive 

GRB HOST DUST EXTINCTION. 







c 2 


C3 


c 4 


X 2 /d.o.f. 


MW 


14.4 


6.52 


2.04 


0.0519 


1.66 


LMC 


4.47 


2.39 


-0.988 


0.0221 


1.19 


SMC 


38.7 


3.83 


6.34 


0. 


1.36 


Linear 


66.2 


4.97 


22.1 


0. 


1.42 


Calzetti 


44.9 


7.56 


61.2 


0. 


1.68 
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Table 2 

Results of fitting to the afterglow SEDs of GRB 00030 1C and GRB 021004 with the Drude approach (see O E]i or various template 

EXTINCTION LAWS. NOTE THAT THE DRUDE APPROACH HAS MORE FREE PARAMETERS THAN THE OTHER APPROACHES. 



Extinction 
Type 


c\ 


C2 


C3 


C4 


Av 
(mag) 


P 


OJy) 




xVMi.o.f. 


GRB000301C 


Drude 


0.025 


0.048 


-2.00 


0. 


0.32 


0.61 


3.99E10 


0.33 


1.98 


MW 










0. 


0.85 


1.04E14 


1.32 


2.64 


SMC 










0.11 


0.62 


4.68E10 


0.64 


1.28 


LMC 










0. 


0.85 


1.04E14 


1.32 


2.64 


Linear 










0.20 


0.51 


1.23E9 


0.47 


0.94 


Calzetti 










0. 


0.85 


1.04E14 


1.32 


2.64 


GRB 021004 


Drude 


0.015 


0.15 


-2.00 


0. 


0.13 


0.78 


6.13E12 


0.47 


1.64 


MW 










0. 


1.06 


8.23E16 


1.53 


2.68 


SMC 










0.15 


0.67 


1.58E11 


0.36 


0.53 


LMC 










0. 


1.06 


8.23E16 


1.53 


2.68 


Linear 










0.26 


0.54 


2.05E9 


0.76 


1.33 


Calzetti 










0.95 


0. 


46.7 


0.80 


1.40 



Table 3 

Results of Drude-fitting to the artificial SED generated by reddening the power-law afterglow F„(/dy) = 5.2 x 10 s (z^/Hz) -0 ' 5 with 

A v = 0.5 MAG EXTINCTION OF MW, SMC, AND CALZETTI-TYPE (SEE FlG.01. 



Reddening ~i c~ 2 c~ 3 cl ~Ay /3 ¥ xV^daia X 1 / N i.o.f. 
Type (mag) (pjy) 



TlW B73 549 2M 0.0514 0301 0499 5.24E8 3.26E-4 4.35E-4 

SMC 39.4 3.89 6.31 0. 0.500 0.501 5.26E8 1.32E-3 1.76E-3 
Calzetti 45.2 7.51 61.7 0. 0.497 0.502 5.17E8 7.98E-4 1.06E-3 
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FIG. 2. — Upper panel: fitting the SED of the afterglow of GRB 000301C (filled black circles) with the SMC (green) and LMC or MW (blue) template extinction 
laws and the "Drude" approach (red; see Eq|2) for the host extinction curve. No extinction is allowed in the MW and LMC models (i.e. the best fit with a MW- or 
LMC-type extinction is given by Ay K, 0): a small amount of Ay would lead to large deviations from the afterglow SED since the 2175 bump prominent in the MW 
and LMC laws is absent in the afterglow SED. Lower panel: comparison of the SMC (green), LMC (blue), MW (Ry =3.1; black) extinction laws with that derived 
from the Drude approach (red). 



8 



i/(GHz) [host frame] 

1E5 5E5 1E6 1.5E6 2E6 2.4E6 




FIG. 3.— Same as Figure|2]but for GRB 021004. 
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FIG. 4. — Upper panel: Drude fits to the observer-frame UBVRUHK "photometry data" artificially-generated by reddening the intrinsic afterglow spectrum 
F„ oc of a burst at z ~ 2 (black line with red crosses superimposed for the observer-frame UBVRUHK bands) with the SMC ("data": cyan squares; Drude fit: 
green line), Calzetti ("data": blue triangles; Drude fit: magenta line), and MW ("data": black circles; Drude fit: red line) extinction laws (with Ay = 0.5 mag for 
each). Lower panel: comparison of the SMC, Calzetti and MW extinction curves (solid lines) with that inferred from the Drude approach (dashed lines). 



